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We consider tlie secular behavior of shell-model configuration (partial) densities. When config- 
uration densities are characterized by their moments, one often finds large third moments, which 
can make suitable parameterization of the secular behavior problematic. We review several param- 
eterizations or models, and consider in depth three specific models: Cornish-Fisher, binomial, and 
modified Breit-Wigner distributions. Of these three the modified Breit-Wigner provides the best 
secular approximation to exact numerical configuration densities computed via full diagonalization 
from realistic interactions. 
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INTRODUCTION AND MOTIVATION 



Neutron-capture rates onto compound nuclear states are usually computed using the statistical Hauser-Feshbach 
formalism^. One important, and uncertain, input into Hauser-Feshbach calculations is the density of excited states, 
t — ■ or level density in addition to the total level density one also would like the exciton (particle- hole) density for 
CN , pre-equilibrium emission 0, |^, |^ . 

The nuclear level density is not trivial to extract experimentally, and there has recently been increasing theoretical 
efforts to compute and characterize the level density; it is beyond the scope of this paper to characterize all approaches, 
although some recent references are 0, 121; llli ilil 112b J>2e JJi J£J • Instead we focus on microscopic models, and 
in particular on the interacting shell model, which accurately describes low-lying spectra and transitions for a broad 
range of nuclides. On the other hand, to extract levels densities from traditional shell model codes requires full 
diagonalization of the Hamiltonian, a computationally forbidding requirement. An alternative to diagonalization 
. is the Monte Carlo path integral technique 0, which is well suited to thermal observables mil El. Although 
reasonable successful, path integral methods are limited to interactions that are free of the "sign problem" |l9ll2Ctl21j . 
Therefore we feel motivated to consider an alternate method based on spectral distribution theory (also known as 
statistical spectroscopy) p3. l2^ . 

Spectral distribution theory computes the moments of the shell-model Hamiltonian. One must then invert the 
moments to find the level density. Note that most other methods to compute the level density also require inversions, 
such as inverse Laplace transform through the saddle-point approximation 11, '12*1 or maximum entropy methods'?]; in 
(~| those approaches, as here, the success of the inversion depends upon the validity of explicit and implicit assumptions 
about the level density. For moment methods one must choose a parameterized "model" for the secular behavior, and 
then adjust the parameters to fit the calculated moments; the choice of model implicitly makes assumptions about 
the higher moments which were not fitted. 



5h One can model the total density of states as a sum of partial or configuration densities [2J, |25| . It is more efficient 

- - 1 and tractable to compute the low moments of many configurations (subspaces) , which one can do directly from the 



two-body matrix elements |23l l26l l27l|. rather than many moments of the entire space. A further advantage is that 
one automatically gets out the particle-hole (exciton) density needed for calculating pre-equilibrium emission. 

The success of a moments-method approach to the level density depends in large part on using a suitable model 
parameterization for the configuration density. Early attempts to find a suitable model for the level density utilized 
Gram-Charlier and Edgeworth series representations in terms of derivatives of an asymptotic density f28|; in essence 
one starts with a Gaussian ,2 2j and expands about it usin g or thogonal polynomials. Fixed- J expansions methods have 
also been developed in terms of Gaussian shapes [23, Hfl • 

Realistic configuration densities often have a large asymmetry, or third moment [2^ l32|: an earlier study demon- 
strated that Gram-Charlier/Edgcworth expansions do poorly for large asymmetries |2j. 

In this paper we first describe the desired features of any distribution used to to model the configuration level 
density. We then review several model parameterizations, focusing in particular on three: Cornish-Fisher, which had 
been shown superior to Gram-Charlier/Edgeworth in |24 |: binomial distributions js^; and finally a new proposal, a 
modified Breit-Wigner distribution. After comparing strengths and weaknesses of these three models, we illustrate 
their performance against exact shell-model calculations, with realistic interactions, of partial densities with a range 
of different asymmetries. We conclude that the modified Breit-Wigner (MBW) has significant advantages over the 
other two distributions. 
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II. CONFIGURATION MOMENTS 



Spectral distribution theory, or nuclear statistical spectroscopy, analyzes many nuclear properties through low- 
lying moments of the Hamiltonian . In this section we review the needed definitions; a somewhat more detail 
discussion is found in 32]. We work in a finite model space M wherein the number of protons and neutrons is fixed. 
If in A4 we represent the Hamiltonian as a matrix H, then all the moments can be written in terms of traces. The 
total dimension of the space is D = trl, and the average is (O) = D~^trO. The first moment, or centroid, of the 
Hamiltonian is E = (H); all other moments are central moments, computed relative to the centroid: 

= ((H - ED , n>2 . (1) 

The width a = y^JjJ^, and one scales the higher moments by the width: 

= ^. (2) 

In addition to the centroid and the width, the next two moments have special names. The scaled third moment 
m*-^^ is the asymmetry, or the skewness; similarly, m^^^-S is the excess (hence a Gaussian has zero excess). 
We use a, /3, 7, . . . to label subspaces. Let 

be the projection operator for the a-th subspace. In this paper we use spherical shell-model configurations to partition 
into subspaces; a configuration is all states of the form, e.g. (131/2)^(151/2)^(0^3/2)^(0^3/2)!^, etc.. One could use 
other group-theoretical partitions but spherical configurations have been the most widely studied for moment methods. 
Now one can introduce partial or configuration densities, 

Pa,{E)^tTP^SiE-li). (4) 

The total density p{E) is just the sum of the partial densities. (Incidentally, one can differentiate between state 
densities, which includes all 2J-I-1 degeneracies in Jz, and level densities, which do not. Be aware, however, that both 
terms are sometimes used interchangeably. While our discussion here is germane to both state and level densities, all 
our specific numerical examples refer to state densities.) 

With projection operators for subspaces Pa in hand, we define partial or configuration moments: the configuration 
dimension is = trP^, the configuration centroid is E^ = D~^tT PaH., while the configuration width CTq and 
configuration asymmetry ma are defined in the obvious ways. An early study of typical configuration moments with 
realistic interactions was done in while more recently we conducted a detailed investigationjs^ . 

One can compute the configuration moments directly from the many-body Hamiltonian, but that is computationally 
unfeasible for large systems. (For our examples herein, however, we do exactly that, generating the many-body 
Hamiltonian matrix via the shell-model code REDSTICK|3J| .) Alternately, one can make use of the available 'analytic' 
formulae for configuration moments 126L 1271 : these still require nontrivial computational effort, especially for third 
and fourth moments, but because they do not require the intermediate step of computing the many-body matrix 
elements, for large systems they are still faster. 

Finally, if one has the a secular density p{E) (which may model either a total density or a configuration density), 
one can compute the moments through integrals rather than traces: 

D = J dEp{E), (5) 

^ = 1 y dEEp{E), (6) 

^(2) ^ ^ 1 dE{E-Efp{E), (7) 

and so on. Formally these integrals are equivalent to the trace definitions. An important question in this paper is 
how accurately the model parameterizations actually reproduce the desired moments, which we will address through 
numerical integration. 
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III. MODELS FOR SECULAR BEHAVIOR 



We want a parameterized model Pmodci(^') for the secular behavior of the (configuration) density of states, with 
the parameters fixed by the low-lying moments of the density. That is, one finds the exact low-lying many-body 
configuration moments, the target moments, from the Hamiltonian and then finds the parameters of the secular model 
that, in principle, will reproduce those target moments. The ideal characteristics for secular behavior are 

1. Non-negative densities. Negative- valued densities are of course unphysical. 

2. Ease of deriving model parameters from target moments. This is important when one has thousands or even 
hundreds of thousands of configuration densities to construct. Furthermore, as discussed below, some 'analytic' 
expressions for moments of model functions are in fact not accurate. 

3. Fixed start- and end-points. While not essential, finite endpoints are useful and reflect the finite range of 
densities in a finite model space. 

Not all models will contain all of these desirable characteristics. We will discuss a total of five models. The first two 
models, Gram-Charlier/Edgeworth and exponential with polynomial argument (EPA), we discuss only briefly. The 
final three- Cornish-Fisher, binomial, and modified Breit-Wigner-we consider in detail. 

Several of these models begin with a Gaussian with centroid E and width a: 



P{E) 



1 



V27rcr2 



exp 



E-E 



(8) 



The Gram-Charlier/Edgeworth distributions |35l | modify a Gaussian by multiplying it by a linear superposition of 
orthogonal Hermite polynomials. This intuitive approach has the highly desireable feature that, given the target 
moments, the parameters are trivial to determine. Unfortunately, these distributions can also have unphysical, non- 
negative densities, which are particularly severe for highly asymmetric distributions. In addition, a previous studvp^l 
showed such distributions simply do poorly in reproducing the secular behavior of realistic densities. 

Another model which generalizes a Gaussian is the exponential with polynomial argument (EPA) distributionjs^ , 
which takes the form p{E) ^ exp{aE + bE^ + cE^ + dE'^) . This distribution is positive definite, but has the problematic 
feature that the moments must be computed numerically. (In practice one uses a look-up table and interpolates (22iJ)- 

For this paper we focus on three models: 

1. Cornish-Fisher. We use the representation found in Ref. "2^, where it was studied and found to be superior to 
Gram-Charlier/Edgeworth distributions. We write it explicitly to show its connection to a Gaussian distribution: 
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FIG. 1: Moment discepancies in Cornish-Fisher (left) and continuous binomial (right) distributions. Each graph corresponds 
to the numerical moments compared with target moments, at different asymmetries mS^\ The plots for Cornish-Fisher are 
with m''*-' = 3.0 only; distributions with higher fourth moments show even larger discrepancies. 



While the Cornish-Fisher distribution is positive definite, and was shown |2J| to be superior to Gram- 
Charlier/Edgeworth distributions, it has some significant drawbacks. Most important is that the relations Hl()|l 
between the Cornish-Fisher parameters and the target moments are approximate^ derived assuming small deviations 
from a Gaussian. 

We checked the consistency of the parameterization (|1U|I by computing moments numerically. That is, given some 
target moments we used relations (I1U|I to obtain the Cornish-Fisher parameters, and then obtained the numerical 
moments using 0. The lefthand side of Fig. shows the discrepancy between target and numerical moments. Not 
only do significant deviations develop beyond |m*-'^^ > 0.5|, the Cornish-Fisher distribution starts to develop unphysical 
bumps at the tails. 

2. Binomial. The binomial distribution [s^ starts from the discrete expansion (1 + A)^ = '^'^ ^ ^ • 1^ 

excitation energies come in discrete steps e, that is. Ex = fee, one can compute the moments for the discrete distribution 
exactly: 



iVeA 




m(4) - 3 = (m(3)) -j-. (11) 
One can obtain a continuous distribution by using gamma functions: 



/pN ^ \-E/e T{Eraaxl(^+ 1) 



(12) 



where i^max — Ne. The distribution can be shifted by adding a simple displacement Smin to the energy. The 
parameters N, A are fitted to the third moment and to either the total dimension d — {1 + A)^ or, through rescaling, 
to the fourth moment. 

In the limit of no asymmetry and N oo the binomial distribution goes to a Gaussian. The binomial distribution 
is positive-definite and has finite endpoints at {E^in,Ne + i?min)- Unfortunately, once again, the relation between 
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TABLE L Comparison of the moments reproduced by modified Breit-Wigner, Cornish-Fisher (C-F) and continuous binomial 
distributions. They are compared against exact calculations, all with GXPFl interaction 38;] on p/-shell ''^Ti for the same 
three configurations shown in Fig. |21 
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the parameters and the target moments, exact for the discrete distribution, is only approximate for the continuous 
distribution. We compared the discrepancy between target moments and numerical moment in the righthand side of 
Fig. n Once again, the approximate formulas (|ll|l is restricted to jm'^'^^l < 0.5. 

3. Modified Breit-Wigner (MBW). Given the problems with the previous models, we abandon generalizations of 
Gaussians and instead propose 

= {E-E,Y,-W^ 

with endpoints given by -Emin , £'max ■ The first four moments can be computed analytically in terms of the four 
parameters; given target moments one can find the parameters by solving four nonlinear algebraic equations (plus 
an overall scale to match the correct dimension); we give the expressions in closed (but nontrivial) form in the 
Appendix. Unlike the Cornish-Fisher or binomial distributions, there is no discrepancy between the target and 
numerical moments. 

The MBW distribution is applicable for 

mt'^) -3 > 1.42 (^m(3))^ -0.52 . (14) 

For third and fourth moments inside the allowed area, one can find real solutions for the MBW parameters. This 
region includes the numerically observed range for realistic interactions and model spaces |33 |. 



A. Direct comparison of model distributions 



To facilitate comparison of the Cornish-Fisher, binomial, and MBW distributions, we make two side-by-side com- 
parisons of these three models along with exact numerical configuration densities; specifically, we computed several 
partial densities for ^''Ti (2 protons and 2 neutrons in a p/-shell valence space) with the realistic interaction GXPFl 
[33 , and chose to use three different configurations with different asymmetries m!^^^ . 

As discussed above, the "analytic" moments of the Cornish-Fisher and continuous binomial distribution are only 
approximate. Table compares "exact" or target moments (shell- model calculations with REDSTICK code 34|) 
and the numerical moments for the Cornish-Fisher and continuous binomial distributions. (The MBW numerical 
moments, as previously discussed, agree with the target moments.) The numerical Cornish-Fisher moments do well 
for centroids and widths, less well for m!^^\ and poorly for m*^^^; the continuous binomial has even larger errors. 

One could in principle adjust the Cornish-Fisher or binomial parameters (with a different, or even an exact, 
parameterization) to get agreement between target and numerical moments, but this would have to be done numerically 
or through a look-up table. Furthermore, the binomial distribution has a much smaller region of applicability than 
the MBW distribution; one cannot find solutions for \m!'^^ \ > 0.5. 
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FIG. 2: Representation of partial densities for three different configurations of Ti in the pf-shell, with GXPFf interaction. 
The modified Breit-Wigner, Cornish-Fisher and Binomial distributions are compared against the exact densities (solid plots 
with binned lines). 



In Fig. 121 we compare the MBW, Cornish-Fisher, and binomial distributions against reahstic configuration densities 
from direct diagonahzation using REDSTICK. Of the three models, the MBW appears to be better. Finally, we are 
able to compute the total level density with the MBW model for '''^Ti in the p/-shell with the GXPFl interaction. 
Fig. O shows such calculation, and the exact calculation with the REDSTICK code, as well. In this example, the 
total level density is a sum of a 100 partial densities. 



IV. SUMMARY 



In order to invert configuration moments to get level densities, one must posit a model distribution whose parameters 
are fixed by the low-lying moments. An important consideration for modeling partial or configuration densities are 
model distributions that perform well for large asymmetries (third moments) which we have previously shown to 
be important [32l|. Continuing previous investigation 24], we have compared several distributions, in particular the 
Cornish-Fisher, continuous binomial, and a modified Breit-Wigner (MBW). We suggest that the MBW, or perhaps 
some variant, may be the most reliable for configuration densities, especially at large asymmetries. 
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FIG. 3: Full '*''Ti level density in the p/-shcll. The modeled level density (sum of MBW partial densities, broken line) is 
compared against the exact calculation (binned line). 
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APPENDIX: MOMENTS EXPRESSIONS FOR THE MODIFIED BREIT-WIGNER DISTRIBUTION 

The modified Breit-Wigner distribution, with dimensions of [Energy]"^, is 

1 {E — EminY {Emax — 



p{E) 



{E - Eof + 



defined on the interval £^min < E < Smax- 
Define the dimensionless integrals 



Io{a,b) 



dx = arctanfe — arctana, 

.t2 + 1 



h{a,b) 



dx. 



Also define the general integral (with an easily proved recursion relation). 



x' + l 



-dx 



so that 
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2k 



Now introduce the convenient dimensionful integral 

rb 



M„= f {E - EoT p{E)dE. 

J a 



(A.1) 



(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 

(A.7) 



We can write this in terms of the dimensionless integrals, using a = (-Bmin — Eo)/W and b = (-Bmax — Eo)/W, 



M„ = W {ln+4 - 2(a + 6)7„+3 + (a^ + 4a6 + b'^)In+2 - 2(a + 6)a6/„+i + aH^In) • 
Finally, the moments in terms of the above expressions are 



dimension D = f p{E)dE = Mq 

J a 



centroid E 



variance a 
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{E- Efp{E)dE = — (Ma - 2A£Mi + AE'^Mq) 
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(Ms - ^AEM2 + SAS^Mi - AE^Mq) 
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(A.IO) 
(A.11) 
(A.12) 



[Mi - AAEMs + 6AE'^M2 - AAE^Mi + AE'^Mq) . (A.13) 
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